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We describe the application of Dyson-Schwinger equations to the calculation of 
hadron observables. The studies at zero temperature (T) and quark chemical po- 
tential (fi) provide a springboard for the extension to finite-(r, /i). Our exemplars 
highlight that much of hadronic physics can be understood as simply a manifes- 
tation of the nonperturbative, momentum-dependent dressing of the elementary 
Schwingcr functions in QCD. 



1 DSE Essentials 

The Dyson-Schwinger equations (DSEs) provide a Poincare invariant, contin- 
uum approach to solving quantum field theories. They are an infinite tower 
of coupled integral equations, with the equation for a particular n-point func- 
tion involving at least one m > rt-point function. A tractable problem is only 
obtained if one truncates the system, and historically this has provided an 
impediment to the application of DSEs: a priori it can be difficult to judge 
whether a particular truncation scheme will yield qualitatively or quantita- 
tively reliable results for the quantity sought. As integral equations, the anal- 
ysis of observables is a numerical problem and hence a critical evaluation of 
truncation schemes often requires access to high-speed computers .0 With such 
tools now commonplace, this evaluation can be pursued fruitfully. 

The development of efficacious truncation schemes in not a purely numer- 
ical task, and neither is it always obviously systematic. For some, this last 
point diminishes the appeal of the approach. However, with growing commu- 
nity involvement and interest, the qualitatively robust results and intuitive 
understanding that the DSEs can provide is becoming clear. Indeed, someone 
familiar with the application of DSEs in the late-70s and early-&Ds might be 
surprised with the progress that has been made. It is now clcarEI that trun- 
cations which preserve the global symmetries of a theory; for example, chiral 
symmetry in QCD, are relatively easy to define and implement and, while it 
is more difficult to preserve. local gauge symmetries, much progress has been 
made with Abelian theoriesB and more is being learnt about non-Abelian ones. 

"The human and computational resources required are still modest compared with those 
consumed in contemporary numerical simulations of lattice-QCD. 
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The simplest truncation scheme for the DSEs is the weak-couphng expan- 
sion. It shows that the DSEs contain perturbation theory; i.e, for any given 
theory the weak-coupUng expansion generates aU the diagrams of perturba- 
tion theory. However, the most important feature of the DSEs is the antithesis 
of this weak-coupUng expansion: the DSEs are intrinsicaUy nonperturbative 
and their solution contains information that is not present in perturbation 
theory. They are ideal for the study of dynamical chiral symmetry breaking 
(DCSB) and confinement in QCD, and of hadronic bound state structure and 
properties. In this application they provide a means of elucidating identifiable 
signatures of the quark-gluon substructure of hadrons. 

Their intrinsically nonperturbative nature also makes them well suited to 
studying QCD at finite-(T, /i), where the characteristics of the phase transition 
to a quark-gluon plasma are a primary subject. The order of the transition, 
the critical exponents, and the response of bound states to changes in these 
intensive variables: all must be elucidated. The latter because there lies the 
signals that will identify the formation of the plasma and hence guide the 
current and future experimental searches. 



1.1 Gluon Propagator 

In Landau gauge the two-point, dressed-gluon Schwinger function, or dressed- 
gluon propagator, has the form 



Qik'^) nli2\ .9 
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i,'l>,.W=(^*,.-^j^,«M:=(j^, (1) 

where n(fc^) is the vacuum polarisation, which contains all the dynamical infor- 
mation about gluon propagation. Studies of the gluon DSE have been reported 
by many authorga with the conclusion that, if the ghost- loop is unimportant, 
then the charge-antiscreening 3-gluon vertex dominates and, relative to the 
free gauge boson propagator, the dressed gluon propagator is significantly en- 
hanced in the vicinity of = 0.l3 The enhancement persists to fc^ ~ 1-2 GeV^, 
where a perturbative analysis becomes quantitatively reliable. In the neigh- 
bourhood of A;^ = the enhancement can be representedO as a regularisation 



''The possibility that Q(k'^) is finite or vanishes at A;^ = is canvassed in these proceedings. 
In the absence of particle-like singularities in the quark-gluonryertex such behaviour is very 
difficult to reconcile with the observable phenomena of QCD.Q A particle-like singularity is 
one of the form {P'^)~°' , a S (0, 1]. In this case one can write a spectral decomposition for 
the vertex in which the spectral densities are non- negative. This is impossible if a > 1. a = 1 
is the ideal case of an isolated, (5-function singularity in the spectral densities and hence an 
isolated, free-particle pole, a S (0, 1) corresponds to an accumulation, at the particle pole, 
of branch points associated with multiparticle production. 
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of 1/fc"' as a_distribution. A dressed-gluon propagator of this type generates 
confinementij and DCSB without fine-tuning. 



1.2 Quark Propagator 

In a covariant gauge the two-point, dressed-quark Schwinger function, or dressed- 
quark propagator, can be written in a number of equivalent forms 

SiP) ■■= ...^ ^ ^ ^.^^ (2) 
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ij ■ p AijP') + B{p'^) 
is the dressed-quark self energy, which satisfies 



i"fpav{p ) + crs{p ) • (3) 



T.{p) = {Z2-l)n-p+ Z4mhn^+ Zi I 

Jq 



^ g^D^,ip-q)^^^S{q)r:{q,p), (4) 



where r"(g;p) is the renormalised dressed-quark-gluon vertex, mbm is the La- 
grangian current-quark bare mass and :— J (i^g/(27r)^ represents mnemon- 
ically a translationally-invariant regularisation of the integral, with A the regu- 
larisation mass-scale. The quark-gluon-vertex and quark wave function renor- 
malisation constants, A^) and Z2(/i^,A^), depend on the renormalisa- 

tion point, /i, and the regularisation mass-scale, as does the mass renormalisa- 
tion constant ^^(/x^, A^) := Z2ifJ,'^ , A^y^ Zi{fM^ , A^). 

The cLuark mass-function is Af(p^) :— B{p'^) / AijP') and, as illustrated 
in Fig. solving the quark DSE using an infrared-enhanced dressed-gluon 
propagator and a dressed-quark-gluon vertex, (p, q) , that does not exhibit 
particle-like singularities at {p — q)^ — 0, one obtains a quark mass-function 
that mirrors the infrared enhancement of the dressed-gluon propagator. These 
results were obtained using the one-loop formula for the running mass, with 
current-quark masses corresponding to 

™lGcV ™lGcV ™lGcV ™lGeV 
^u/d ™s ™c ™6 

6.6 MeV 140 MeV 1.0 GeV 3.4 GeV. ^' 

The quark DSE was also solved in the chiral limit, which in QCD is ob- 
tained by setting the Lagrangian current-quark bare mass to zeroQ One ob- 
serves immediately that the mass-function is nonzero even in this case. That 



'^One aspect of confinement is the absence of quark and gluon production thresholds in 
colour-singlet-to-singlet <S-matrix amplitudes. This is ensured if the dressed-quark and - 
gluon propagators do not have a Lehmann representation. 



3 



3 



> 



b-quark 
c-quark 
s-quark 
u/d-quark 
chiral limit 
M(p) = p 




2 3 4 5 

p (GeV) 

Figure 1: Dressed-quark mass-function obtained in solving the quark DSE. 



is DCSB: a momentum-dependent quark mass, generated dynamically, in the 
absence of any term in the action that breaks chiral symmetry explicitly. This 
entails a nonzero value for the quark condensate in the chiral limit. That 
M(p^) 7^ in the chiral limit is independent of the details of the infrared- 
enhanced dressed-gluon propagator. 

Figure |l| illustrates that for light quarks {u, d and s) there are two distinct 
domains: perturbative and nonperturbative. In the perturbative domain the 
magnitude of M{p'^) is governed by the the current-quark mass. For < 
1 GeV^ the mass-function rises sharply. This is the nonperturbative domain 
where the magnitude of M(p^) is determined by the DCSB mechanism; i.e., the 
enhancement in the dressed-gluon propagator. This emphasises that DCSB is 
more than just a nonzero value of the quark condensate in the chiral limit! 

The solution of — Af^(p^) defines a Euclidean constituent-quark mass, 
M^. For a given quark flavour, the ratio /ra^^ is a single, quantitative 
measure of the importance of the DCSB mechanism in modifying the quark's 
propagation characteristics. As illustrated in Eq. (|^), 
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(6) 



this ratio provides for a natural classification of quarks as either light or heavy. 
For light-quarks the ratio is characteristically 10-100 while for heavy-quarks 
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it is only 1-2. The values of this ratio signal the existence of a characteristic 
DCSB mass-scale: M^. At > the propagation characteristics of a flavour 
with < are altered significantly by the DCSB mechanism, while for 
flavours with m'j 3> it is irrelevant, and explicit chiral symmetry breaking 
dominates, ft is apparent from the figure that ~ 0.2GeV~ Aqcd- This 
forms the basis for simplifications in the study of heavy-meson observablesQ 

1.3 Hadrons: Bound States 

The properties of hadrons can be understood in terms of their substructure by 
studying covariant bound state equations: the Bethe-Salpeter equation (BSE) 
for mesons and the covariant Fadde'ev equation for baryons. The mesons have 
been studied most extensively and their internal structure is described by a 
Bethe-Salpeter amplitude obtained as a solution of 

[^H{k;P)],^ = r [xH{q:P)]srKZ{q.k-P), (7) 

Jq 

where XH{q\P) := S{q+)V H[q\ P)S{q-); = diag(5„(g), ^d((z), 5,((z), . . .); 

= q + r]p P, g_ — q — (1 — rjp) P, with P the total momentum of the bound 
state and observables independent of rjp] and r,. . . ,m represent colour-, Dirac- 
and flavour-matrix indices. The amplitude for a pseudoscalar bound state has 
the form 

iEH{k;P)+j-PFH(k;P) (8) 

+ j-kk-PGHik;P) + (J^uk^P,HHik-P) , 

where is a flavour matrix that determines the mesonic channel under con- 
sideration; e.g., T^'^ := (1/2) (A^ + iA^), with {A^ j = 1 ... 8} the Gell-Mann 
matrices. In Eq. (^, K is the renormalised, fully-amputated, quark-antiquark 
scattering kernel and important in the successful application of DSEs is that 
it has a systematic skeleton expansion in terms of the elementary, dressed- 
particle Schwinger functions; e.g., the dressed-quark and -gluon propagators. 
This particular expansionfj with its analogue for the kernel in the quark DSE, 
provides a means of constructing a kernel that, order- by-order in the number 
of vertices, ensures the preservation of vector and axial- vector Ward-Takahashi 
identities. 

In any study of meson properties, one chooses a truncation for K. The 
BSE is then fully specified and straightforward to solve, yielding the bound 



TH{k:P)^T"^5 
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Figure 2: First two orders in a systematic expansion of the quark-antiquark scattering kernel. 
In this expansion, the propagators are dressed but the vertices are bare. 

state mass and amplitude. The "ladder" truncation of K combined with the 
"rainbow" truncation of the quark DSE [F^ 7^ in Eq. (|^)] is the sinmlest 
and most often used. The expansion of Fig. 2 provides the explanatioriEl for 
why this Ward-Takahashi identity preserving truncation is accurate for flavour- 
nonsinglet pseudoscalar and vector mesons: there are cancellations between the 
higher-order diagrams. It also shows why it provides a poor approximation in 
the study of scalar mesons, where the higher-order terms do not cancel, and 
for flavour-singlet mesons, where it omits timelike gluon exchange diagrams. 



2 A Mass Formula 

The dressed-axial- vector vertex satisfies a DSE who se kernel is,-^, and because 
of the systematic expansion described in Sec. 1.3 it foUowsQ that the axial- 
vector Ward-Takahashi identity (AV-WTI): 



5~'(fc+)75^ + 75^5-i(fc_) _ M(^)Ff (fc;P) -Ff (fc;P)M(^) , 



(9) 



[M(^) = diag(TO^, m^', m^, . . .) is the current-quark mass matrix] is satisfied in 
any thoughtful truncation of the DSEs. That entails many important resultsQ 



1) The axial- vector vertex has a pole at 
the leptonic decay constant: 



ijj whose residue is Jh, 



IhP^. = Z2 



-tr 



T") j5l^Siq+)THiq;P)Siq-) 



(10) 
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with the trace over colour, Dkac and flavour indices. 
2) In the chiral limit 



fHEnik; 0) = BoiP) , FR{k; 0) + 2 ///^^(fc; 0) = Aoi^) , 

Gnik; 0) + 2 fHCnik; 0) = 2AJ,(fc2) , i7;j(fc; 0) + 2 /Hi?H(fc; 0) = , 

(11) 

where ylo(fc^) and BQ{k^) are the solutions of Eq. in the chiral limit, and 
Fr, Gr and Hr are calculable functions in r|^. This shows that when chiral 
symmetry is dynamically broken: 1) the flavour- nonsinglet, pseudoscalar BSE 
has a massless solution; 2) the Bethe-Salpeter amplitude for the massless bound 
state has a term proportional to 75 alone, with the momentum-dependence 
of i?/f(fc;0) completely determined by that of Bo{k^), in addition to terms 
proportional to other pseudoscalar Dirac structures that are nonzero in general; 
and 3) the axial- vector vertex, r|^(/c;P), is dominated by the pseudoscalar 
bound state pole for ~ 0. The converse is also true. 

3) The pseudoscalar vertex also has a pole at — —niljj whose residue is 

irH = Zij itr [(T«) '755(9+ (g; P)5(g_ )J . (12) 

4) There is an identity between the residues of the pseudoscalar meson 
pole in the axial-vector and pseudoscalar vertices that is satisfied independent 
of the magnitude of the current quark mass: 

/h ml = rHMH, Mh := traavour [M(^) {t^, (T«)'}] . (13) 
2. 1 Corollaries 

Equation (|l3|) is a mass formula for flavour-octet pseudoscalar mesons. For 
small current-quark masses, using Eqs. (|) and (|ll|), Eq. (|l|) yields 

r% = -^ , - {m)l Z^{^i\K^) t troirac [5™=o(<z)] , (14) 

Ih Jq 

where (^g)" is the chiral limit vacuum quark condensate, which is renormalisation- 
point dependent but independent of the gauge parameter and the regularisa- 
tion mass-scale:El {<iqYl=,i gcV ~ ^ (0.241 GeV)^ . Now one obtains immediately 
from Eqs. (13) and (|lj) 

f^ml = - K + + O (ml) (15) 

fK^^K^ = - K + + O {rhl) , (16) 
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which exemphfy what is commonly known as the Gell-Mann-Oakes-Renner 
relation, [rhq is the renormalisation-point-independent current-quark mass.] 
In a typical, model calculationD Eq. ( |l5| ) is accurate to 7%, while Eq. (16) 
receives corrections of 43% from O (to^) and higher, see Fig. 3. 

In the heavy-quark limit, ?7p = 1 in Eq. (|^) and the heavy- meson velocity 
(Vfj,) and binding energy (E) are defined via: P := m^fz;^, = —1, Mh '■= 

rE 



(Mq + E), with Mh the heavy-meson mass and Mq w Mq w mg; see Sec. 1.2 
In this caseQ, at leading order in 1/Mh, 

Siq + P)-'-'-^, (17) 
2 q ■ V — E 

THiq;P) = /M^rH{q;P), (18) 

where the canonical normalisation condition for Tniq'jP) is independent of 
Mh- Using Eqs. (p^) and (pi) in Eq. (fol) one obtains 



fn cx (19) 
and this along with Eqs. ([l^ and (|l^) yields 

Af// oc rhg . (20) 
A model studyi shows Eq. (pOl) to be valid for rhq > rhg, and this is confirmed 



by data. Fig. 3. However, both calculations and available data suggest that 
Eq. ( |l9| ) is not manifest until rhq > The- 



3 Finite T and /i 

The study of QCD at finite temperature and baryon number density proceeds 
via the introduction of the intensive variables: temperature and quark chemical 
potential. These are additional mass-scales, with which the coupling can run 
and hence, for T » Aqcd and/or ^ » Aqcd, Q;s(Q^ — 0, T, fi) ^ 0. It follows 
that at finite temperature and/or baryon number density there is a phase 
of QCD in which quarks and gluons are weakly interacting, irrespective of the 
momentum transfer; i.e., a quark-gluon plasma. Such a phase of matter existed 
approximately one microsecond after the big-bang. In this phase confinement 
and DCSB are absent and the nature of the strong interaction spectrum is 
qualitatively different. 

Nonperturbative methods are necessary to study the phase transition, 
which is characterised by qualitative changes in order parameters such as the 
quark condensate. One widely used approach is the numerical simulation of 
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Fignre 3: Pseudoscalar meson mass as a function of the mass of the heaviest constituent, 
rhqQ Only the tt does not Uo on the same straight line. 



finite temperature lattice-QCD, jjiith the first simulations in the early eighties 
and extensive efforts since thenE£l The commonly used quenched approxima- 
tion is inadequate for studying the phase diagram of finite temperature QCD 
because the details of the transition depend sensitively on the number of active 
(light) flavours. It is therefore necessary to include the fermion determinant. 

That is even more important in the presence of /i, which modifies the 
fermion piece of the Euclidean action: "/■d + m—^j-d — "fiH + m, and thus 
the fermion determinant acquires an explicit imaginary part. The fi QCD 
action being complex entails that the study of finite density is significantly 
more difflcult than that of finite temperature. Simulations that ignore the 
fermion determinant at fi ^ encounter a forbidden region, which begins at 
H = m7r/2p3 and since m.^. — > in the chiral limit this is a serious limitation, 
preventing a reliable study of chiral symmetry restoration. J'he phase of the 
fermion determinant is essential in eliminating this artefactLiJ 

The contemporary application of DSEs at finite temperature and chemical 
potential is a straightforward extension of the T = = /i studies. The direct 
approach is to develop a finite-T extension of Ansdtze for the dressed-gluon 
propagator. The quark DSE can then be solved and, having the dressed-quark 
and -gluon propagators, the response of bound states to increases in T and fi 
can be studied. As a nonperturbative approach that allows the simultaneous 
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study of DCSB and confinement, the DSEs have a significant overlap with lat- 
tice simulations: each quantity that can be estimated using lattice simulations 
can also be calculated using the DSEs. This means they can be used to check 
the lattice simulations, and importantly, that lattice simulations can be used 
to constrain their model-dependent aspects. Once agreement is obtained on 
the common domain, the DSEs can be used to explore phenomena presently 
inaccessible to lattice simulations. 



3.1 Fimte-{T,n) Quark DSE 

The renormalised dressed-quark propagator at finite- (T,/z) has the form 

S(p,u}k) = (21) 

■ pA{p,Ulk) +l^4,UlkC{p,LVk) +B{p,UJk) 

= -ij -pcTAip, Cbk) - i^ALOkCFciv, Wfe) + fTB(p, Wfe) (22) 

where ujk '■= tOk+in with tuk = (2fc-|-l) ttT, the fermion Matsubara frequencies, 
A; € Z. The complex scalar functions: A{p,u)k), B{p,LJk) and C{p,u>k) satisfy: 
J^{p, tjfe)* = J-{p, cj-jt-i) , = A, B,C, and although not explicitly indicated 
they are functions only of and d)^. 
S{p,u)k) satisfies the DSE 

S~^{p, Qk) = Z2 n ■ p + Z2 (174 ojk + Wbm) + S'(p, tDfc) ; (23) 

where the regularised self energy is 

s'(p, uik) = n • p'^'a{P' ^k) + m ^'c{p, ^k) + Sb(p, cbk) , (24) 

4 1 
^'j^{p,^^k)= / -5^-D,,i.(p- g,Wfe - wz)-tr [P:r7;,5'(g,(:;i)ri,(g,W(;p,(I;fe)] , 

Jl,q * 

(25) 

C ■■= ^ (C?^ and Va := -{Z^/P')il ■ P, Vb ■= ^i, Vc := 

-(■^l/Wfc)i74. 

The finite- (T, /i). Landau-gauge dressed-gluon propagator has the form 

g^D^, {p, n) = P^, (p, n) Af if, O) + Pj, ip) Ag (p, fi) , (26) 

{0; n and/or i/ = 4, 

S,-'^; ,,^^ = ^,J = W^ (27) 

with Pff^ip) + P/^^{p,P4) = Sf„y - Pp.PvIYL=y PaPa, H,V = I, ■ ■ ■ A- 
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In studying confinement one cannot assume that the analytic structure of 
a dressed propagator is the same as that of the free particle propagator: it must 
be determined dynamically. The pk '■= (p, d)fc)-dependence of A and C is quali- 
tatively important since it can conspirepsrith that of B to eliminate free-particle 
poles in the dressed-quark propagatorEj In this case the propagator does not 
have a Lehmann representation so that, in general, the Matsubara sum cannot 
be evaluated analytically. More importantly, it either complicates or precludes 
a real-time formulation of the finite temperature theory, which makes the study 
of nonequilibrium thermodynamics a very challenging problem. In addition, 
the pfc-dependence of A and C can be a crucial factor in determining the be- 
haviour of bulk thermodynamic quantities such as the pressure and entropy, 
being responsible for these quantities reaching their respective ultrarelativistic 
limits only for very large values of T and fi. It is therefore important in any 
DSE study to retain A{pk) and C{pk), and their dependence onpk- 

3.2 Phase Transitions and Order Parameters 

One order parameter for the chiral symmetry restoration transition is well 
known - it is the quark condensate, defined via the renormalised dressed-quark 
propagator, Eq. (Q). An equivalent order parameter is 

X ■.= ReBo{p^Q,Cjo), (28) 

which makes it clear that the zeroth Matsubara mode determines the charac- 
ter of the chiral phase transition. An order parameter for confinement, valid 
for both light- and heavy-quarks, was introduced in Ref. |l^. It is a single, 
quantitative measure of whether or lupt a Schwinger function has a Lehmann 
representation, and it has been usedo to striking effect in QED3. 

3.3 Study at {T ^0,^1 = 0) 

Deconfinement and chiral symmetry restoration have been studiedEl in a DSE- 
model of two-ligkt-flavour QCD. The quark DSE was solved using a one- 
parameter modelEj dressed-gluon propagator, which provides a good descrip- 
tion of TT and /9-meson observables at T = = /i. The transitions are coincident 
and second-order at a critical temperature of w 150 MeV, with the same 
estimated critical exponents: (3 = 0.33 ± 0.03. Both the pion mass, m^r, and 
the pion leptonic decay constant, /^r, are insensitive to T for T < 0.7 Tc. How- 
ever, as T — > Tc, the pion mass increases substantially, as thermal fluctuations 
overwhelm quark-antiquark attraction in the pseudoscalar channel, until, at 
T = Tc, /tt ^ and there is no bound state. These resultSj-confirm those of 
contemporary numerical simulations of finite-T lattice-QCDta. 
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3-4 Complementary study at (T = 0, /i 7^ 0) 

The behaviour of this model at /i ^ has also been explored0 The model 
employs a dressed-gluon propagator with the simple form 



fc2 



16 2 

TT 

9 



1 



-[feV(4m?)] 



fc2 



(29) 



whereB rrit = 0.69 GeV is a mass-scale that marks the boundary between the 
perturbative and nonperturbative domains, and the quark DSE was solved in 
rainbow approximation. The solution has the form 



(30) 



where := (p, W[^]), with wj^j :— + i/i. There are two distinct types of 
solution: a Nambu-Goldstone mode characterised by gbq ^ 0; and a Wigner- 
Weyl mode characterised hy cjbo = 0- 

The possibility of a phase transition between the two modes is explored 
by calculating the relative stability of the different phases, which is measured 
by the difference between the tree-level auxiliary-field effective-action: 



2N;Nc 



(31) 



fl 









where A and C represent the solution of Eq. ( p3| ) obtained when = 0; 
i.e., when DCSB is absent—. This solution exists for all ji. B{ii) defines a 
dependent "bag constant" E3 It is positive when the Nambu-Goldstone phase 
is dynamically favoured; i.e., has the highest pressure, and becomes negative 
when the Wigner pressure becomes larger. Hence the critical chemical potential 
is the zero of which is /ic = 375 MeV. The abrupt switch from the Nambu- 
Goldstone to the Wigner- Weyl mode signals a first order transition. 

The chiral order parameter increases with increasing chemical potential 
up to /.tc, with x(/^c)/x(0) ~ 1.2, whereas the deconfincmcnt order parameter, 
k(/x), is insensitive to increasing /i. At /ic they both drop immediately and dis- 
continuously to zero, as expected of coincident, first-order phase transitions. 
The increase of X with /i is a necessary consequence of the momentum depen- 
dence of the scalar piece of the quark self energy, i3(p[^])tj The vacuum quark 
condensate behaves in qualitatively the same manner as %. 

Even though the chiral order parameter increases with yu, m,r decreases 
slowly as jjL increases, with TO^(/Lt « 0.7/Zc)/to^(0) w 0.94. At this point 
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begins to increase although, for ii < iic, tii-^[^) does not exceed m7r(0), 
which precludes pion condensation. The behaviour of results from mutually 
compensating increases in and {m''j^{qq) (^) tt ■ /tt is insensitive to the chemical 
potential until ^ ~ 0-7 Hc, when it increases sharply so that /tt (Atf!^ ) //tt (/^ ~ 
0) « 1.25. The relative insensitivity of m^r and /tt to changesin fi, until very 
near ^c, mirrors the behaviour of these observables at finite-TEj For example, 
it leads only to a 14% increase in the n ^ decay width at « 0.9 ^c- The 
universal scaling conjecture of Ref. |^ is inconsistent with the anticorrelation 
observed between the /i-dependence of and m^. 

Comparing the /i-dependence of /^r and m^r with their T-dependence, one 
observes an anticorrelation; e.g., at /x = 0, falls continuously to zero as T 
is increased towards Tc ~ 150 MeV. This is a necessary consequence of the 
momentum-dependence of the quark self energy. 

Note A: In calculating these observables one obtains expressions for or 
and hence the natural dimension is mass-squared. Therefore their 
behaviour at finite T and /i is determined by 

Re(.;f^])~[^2y2„^2]^ (32) 

where the T-dependence arises from the introduction of the fermion Mat- 
subara frequency: p4 — > (2fc + 1)'kT. Hence when such a quantity de- 
creases with T it will increase with /i, and vice-versacil 

The confined-quark vacuum consists of quark-antiquark pairs correlated 
in a scalar condensate. 

Note B: Increasing /i increases the scalar density: {—{qq)). This result is an 
expected consequence of confinement, which entails that each additional 
quark must be locally paired with an antiquark thereby increasing the 
density of condensate pairs as /i is increased. 

For this reason, as long as /z < /ic, there is no excess of particles over anJtipar- 
ticles in the vacuum and hence the baryon number density remains zerop i.e., 
p^'^ = , /i < tJLc- This is just the statement that quark-antiquark pairs 
confined in the condensate do not contribute to the baryon number density. 

The quark pressure, P"+'^[/i], can be calculatecO and one finds that af- 
ter dcconfinement it increases rapidly, as the condensate "breaks-up", and 
an excess of quarks over antiquarks develops. The baryon-number density, 
pV''^ = (l/3)aP"+''/a/i, also increases rapidly, with 

p^+'^(M«2Ai,) ^3po, (33) 

where po = 0-16 fm^^ is the equilibrium density of nuclear matter. For compaiij 
ison, the central core density expected in a 1.4 neutron star is 3.6-4.1 po" 
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Finally, ed fi ^ 5/ic, the quark pressure saturates the ultrarelativistic limit: 
pu+d ^ ^4/(-27r2), and there is a simple relation between baryon-density and 
chemical-potential: 

Va.>5mc, (34) 

3 TT 

so that p^^^'^^ (5/ic) 350/30- Thus the quark pressure in the deconfined 
domain overwhelms any finite, additive contribution of hadrons to the equation 
of state. That was anticipated in Ref. where the hadron contribution was 
neglected. This discussion suggests that a quark-gluon plasma may be present 
in the core of dense neutron stars. 



3.5 Simultaneous study of (T =^ 0, fi ^ 0) 

This is a difficult problem and the most complete study@ to date employs 
a simple Ansatz for the dressed-giuon propagator that exhibits the infrared 
enhancement suggested by Ref. ||] : 



g'D,M ^k) - [6,. - j^f^J Stt^ ^ 4o S'ip) , (35) 

with Qk = SfcTrT, the boson Matsubara frequency. As an infrared-dominant 
model that does not represent well the behaviour of Df^^{p, fi^) away from 

« 0, some model-dependent artefacts arise. However, there is significant 
merit in its simplicity and, since the artefacts are easily identified, the model 
remains useful as a means of elucidating many of the qualitative features of 
more sophisticated Ansdtze. _ 

With this model, using the rainbow approximation, the quark DSE iaa 

S~^{p,ujk) = S^'^{p,LUk) + 2'ri^-iuS{p,ujk)lu ■ (36) 

A simplicity inherent in Eq. ( |3^ ) is now apparent: it allows the reduction of 
an integral equation to an algebraic equation, in whose solution many of the 
qualitative features of more sophisticated models are manifest. 

In the chiral limit Eq. ( ^6| ) reduces to a quadratic equation for B{pk), which 
has two qualitatively distinct solutions. The Nambu-Goldstone solution, with 



B{pk) = <; V '/ - -i^fc ' ^<pI) < - (37) 



otherwise 




,2 



Re(p2) < !L 
, otherwise , 
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describes a phase of this model in which: 1) chiral symmetry is dynamicaUy 
broken, because one has a nonzero quark mass-function, B{pk), in the absence 
of a current-quark mass; and 2) the dressed-quarks are confined, because the 
propagator described by these functions does not have a Lehmann representa- 
tion. The alternative Wigner solution, for which 



S(pfc) ^ , C{pk) = i (^1 + ^1 + 1^ j , (39) 

describes a phase of the model with neither DCSB nor confinement. 

The relative stability of the different phases is measured by a (T, /i)- 
dependent vacuum pressure difference, which in the chiral limit is 

6(T,/i)= (40) 
r?4 2N,Nf^ ' dyy' jpe {2pf) - Re (^-^) - In \plCipi) 

with: T = T/rj, p, — Imax is the largest value of I for which wf^^^ < 

(1/4) + ^^ and this also specifies w/^^^ , = —Cuf.pi = {y,u!i+ip). The 
condition B(T, /i) = defines the phase boundary in the (/z, T)-plane. 

Again, the deconfinement and chiral symmetry restoration transitions are 
coincident. For /i = the transition is second order and the critical tempera- 
ture is r° = 0.159 77, which using the value of 77 = 1.06 GeV obtained by fitting 
the TT and p masses corresponds to — 0.170 GeV. This is only 12% larger 



than the value reported in Sec. 3.2, and the order of the transition is the same. 



For any /i 7^ the transition is first-order. For T — the critical chemical 



potential is /ij? = 0.3 GeV, which is « 30% smaller than the result in Sec. 3.4 
The discontinuity in the order parameters vanishes as fi 0. 

The quark pressure, Pq, is calculated easily in this model. Confinement 
Pg = in tl 

approaches 



means that in the confined domain. In the deconfined domain it 



Pur := N^Nf^ (^p^ + 2n^p^f^ + ^^r^T*) , (41) 

the ultrarelativistic, free particle limit, at large values of T and p. The ap- 
proach to this limit is slow, however. For example, at T 0.3 ~ 2T'", or 
p ~ 1.0 ~ S/i^, Pq is only 0.5 Pur- A qualitatively similar result is observed in 
numerical simulations of finite-T lattice-QCD^J This feature results from the 
persistence of the momentum dependent modifications of the quark propagator 
into the deconfined domain, and predicts that there is a "mirroring" of finite-T 
behaviour in the dependence of the bulk thermodynamic quantities on /i. 
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3.6 TT and p properties 

The model discussed in the last section has been used0 to study the (T, /i)- 
dependence of tt and p properties, and to elucidate other features of models 
that employ a more sophisticated ^nsate for the dressed-gluon propagator. For 
example, the vacuum quark condensate takes the simple form 



- (99) = £' dy Re ^ I - y^ - Cuf^ : (42) 

at r = = /i, {-{qq)) = 'n^/iSOir'^) = (O.llTy)^ {-{qq)) decreases with T but 
increases with p, up to a critical value of Pc(T) when it drops discontinuously to 



zero, in agreement with the behaviour reported in Sees. |3.3| and 3^, see Note B. 
This vacuum rearrangement is manifest in the behaviour of the necessarily- 
momentum-dependent scalar part of the quark self energy, B{pk). 

The leptonic decay constant also has a simple form in the chiral limit: 

.n = V^'-^f 'g I (l + 4 - 4^f - ^ A?) . (43) 

Characteristic in Eqs. ( ^ ) and ( ^ ) is the combination p^ — ujf, see Note A. 
Without calculation, Eq. (^) indicates that will decrease with T and in- 
crease with p. This provides a simple elucidation of the results described in 



Sees. 3.3 and 3.4 



The (T, /^)-response of the tt and p masses is determined by the BSE 
VuiPk-, Pi) = Re {tm S{p., + \Pi) TM{p^■, Pi) S{p, - \Pi) 7^} , (44) 

where P( :— (P, ^e), with the bound state mass obtained by considering Pi=o. 
The TT equation admits the solution 



r.(^o) =75(«^i+7-^^2j (45) 

and the calculated (T, /i)-dependence of the mass is depicted in Fig. 4. 

For the p-meson there are two components: one longitudinal and one trans- 
verse to P. The solution of the BSE has the form 



(46) 
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Figure 4: Mp+ and function of fj, for T = 0,0.1. On the scale of this figure, 

is insensitive to this variation of T. The current-quark mass is m = 0.011 »y, which for 

T] = 1.06 GeV yields Mp+ = 770 MeV and ni,r = 140 MeV at T = = ^t. 



where Op^ labels the longitudinal and 9p_ the transverse solution. The eigen- 
value equation obtained from (^J) for the bound state mass, Mp±, is 

2 

^ Re {as{ujl+ - '-M^p^f - [±4+ - \m'p±\ ^v{ojI+ - \M'p±f} = 1 . 

(47) 

The equation for the transverse component is obtained with [— o^q^ — 
(1/4)M2_] in @. Using the chiral-hmit solutions, Eq. @, one obtains 
immediately that 

Mp_ = 1 independent of T and fi. (48) 

Even for nonzero current-quark mass, Mp_ changes by less than 1% as T and 
fi are increased from zero toward their critical values. Its insensitivity is con- 
sistent with the absence of a constant mass-shift in the transverse polarisation 
tensor for a gauge-boson. 

For the longitudinal component one obtains in the chiral limit: 



= ' ~~ ^'T') ■ (49) 



The characteristic combination [/i^ — tt^T^] again indicates the anticorrelation 
between the response of Mp+ to T and its response to fi, and, like a gauge- 
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boson Debye mass, that M^^ rises linearly with for fi = 0. The m ^ 
solution of Eq. (^7|) for the longitudinal component is plotted in Fig. 4: Mp+ 
increases with increasing T and decreases as /i increases. 

Equation (|7|) can also be applied to the 0-meson. The transverse com- 
ponent is insensitive to T and /z, and the behaviour of the longitudinal mass, 
M^+, is qualitatively the same as that of the p- meson: it increases with T and 
decreases with /x. Using r/ = 1.06 GeV, the model yields M^± — 1.02 GeV for 
nis = 180 MeV at T = = ^i. 

In a 2-flavour, free-quark gas at T = the baryon number density is 
Pb = 2fi'^ / {3tt^) , by which gauge nuclear matter density, po = 0.16 fm^"^, 
corresponds to p — po :— 260 MeV = 0.245^7. At this chemical potential the 
algebraic model yields 

Mp+{po) » 0.75Mp+{p = 0) , M^+{po) « OMM^+{p = 0) . (50) 



The study summarised in Sec. 3.4 indicates that a better representation of the 
ultraviolet behaviour of the drcssed-giuon propagator expands the horizontal 
scale in Fig. 4, with the critical chemical potential increased by 25%. This 
suggests that a more realistic estimate is obtained by evaluating the mass at 
p'q — 0.20 77, which yields 

Mp+ip'o) » 0.85Mp+{p = 0) , M^+{p'o) « 0.90M^+(^ = 0) ; (51) 



a small, quantitative modification. The difference between Eqs. ( |50| ) and (51) 
is a measure of the theoretical uncertainty in the estimates in each case. At 
the critical chemical potential for T = 0, A/p+ « 0.65 Mp+{p = 0) and Af0+ « 
OmM^+{p = 0). 

This simple model preserves the momentum-dependence of gluon and quark 
dressing, which is an important qualitative feature of more sophisticated stud- 
ies. Its simplicity means that many of the consequences of that dressing can be 
demonstrated algebraically. For example, it elucidates the origin of an anticor- 
relation, found for a range of quantities, between their response to increasing 
T and that to increasing p. And the (T, /i)-dependence of {—{qq)) and /,rj un- 
derstood algebraically, is opposite to that observed for rrtp^, hence the scaling 
law conjectured in Rcf. |po| is inconsistent with this calculation, as it is with 
others of this type. 



4 Concluding Remarks 

This contribution illustrates the contemporary application of Dyson- Schwinger 
equations to the analysis of observable strong interaction phenomena, high- 
lighting positive aspects and successes. Many recent, interesting studies have 
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been neglected: caJrulations of the cross section for diffractive— vector meson 
electroproductionl^J the electric-jdipole moment of the /o-mesonpiand the elec- 
tromagnetic pion form factor;c3 an exploration of 77-77' mixing;L£l and others 
reviewed in Ref. ||2^ ]. However, a simple enquiry of 

" http: / /xxxdanl.gov/find/hep-ph ' 
with the keywords: "Dyson-Schwinger" or "Schwingcr-Dyson" , will provide a 
guide to other current research. 

In all phenomenological applications, modelling is involved, in particular 
of the behaviour of the dressed Schwinger functions in the infrared. (The ultra- 
violet behaviour is fixed because of the connection with perturbation theory.) 
This is tied to the need to make truncations in order to define a tractable 
problem. Questions will always be asked regarding the fidelity of this mod- 
elling. The answers can only come slowly as, for example, more is learnt about 
the constraints that Ward Identities and Slavnov- Taylor identities in the. the- 
ory can provide. That approach has been particularly fruitful in QEDp and 
already in the development of a systematic truncation procedure for the ker- 
nel of the quark DSE and meson BSEEj In the meantime, and as is common, 
phenomenological applications provide a key to understanding which elements 
of the approach need improvement: the approach itself must also be explored 
under extreme conditions. 
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